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ABSTRACT 

We present an X-ray stacking analysis of a sample of 38 submillimeter galaxies with (z) =2.6 
discovered at > 4cr significance in the Lockman Hole North with the MAMBO array. Wc find 
a 5(7 detection in the stacked soft band (0.5-2.0 keV) image, and no significant detection in the 
hard band (2.0-8 keV). We also perform rest-frame spectral stacking based on spectroscopic and 
photometric redshifts and find a ~ 4cr detection of Fe Ka emission with an equivalent width 
of EW > IkeV. The centroid of the FeKa emission lies near 6.7 keV, indicating a possible 
contribution from highly ionized FeXXV or FeXXVI; there is also a slight indication that the 
line emission is more spatially extended than the X-ray continuum. This is the first X-ray analysis 
of a complete, flux-limited sample of SMGs with statistically robust radio counterparts. 

Subject headings: Galaxies: Submillimeter, Galaxies: active. Galaxies: formation. X-rays: galaxies 



1. Introduction 

Submillimeter galaxies (SMGs) are distant star- 
forming systems with tremendous infrared lumi- 
nosities (LiR [8-1000 ^m] > 10^2 Lq). In the 
(sub)millimeter waveband they are observable out 
to high redshifts due to the strong negative K- 
correction in the Raylcigh- Je ans regime of thei r 



thermal spectrum (see, e.g., Blain et al] 20021) 



The p revalence of SMGs at z > 1 (jChapman et al 



20051 ) in combination with their high rates of 
dust-obscured star-formation imply that they may 
be responsible for the production of a signifi- 
cant fraction of all the s tellar mass in p r esent 
day g alaxies. X-ray (I Alexander et al. 
2005b) and mid-infrared (Valiantc et al 



2003, 



2007 



Menendez-Delmestre et al.ii2007. . . 2009; .Pope et al 
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2008f ) spectroscopy shows that SMGs frequently 
contain active galactic nuclei (AGN) as well as 
powerful starbursts. This connection between star 
formation and accretion at high redshift may help 
explain the black hole mass-bul ge mass relation 



in pre sent-day galaxies (e.g., [Alexander et al 
2005al ). However, it remains hard to determine 



the relative importance of accretion and star for- 
mation for the SMG population as a whole because 
of the challenge of assembling large, statistically 
unbiased SMG samples. 

Studying the X-ray properties of SMGs is diffi- 
cult for two main reasons. First, the X-ray coun- 
terparts to SMGs are extremely faint. The count 
rate is so low that even the deepest Chandra 
and XMM — Newton spectra of SMGs cannot 
resolve features that serve as sensitive diagnos- 
tics of the physical conditions inside galaxies, 
like the FeKa emission line. FeKa emission is a 
ubiquitous feature in spectra of optically-selected 



AGN up to z c± 3 ( e.g.. iBrusa et al. 1 12005 



Chaudharv et al. 2010l : llwasawa et al.l 2011 ) , but 



the FeKa emission properties of SMGs remain 
consid erably more uncertain (lAlexander et al. 
2005bl) . Second, the requirement that SMGs need 



radio, (sub)millimcter, or mid-IR counterparts 
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capable of nailing down their positions in high- 
resolution X-ray maps can lead to concessions of 
inhomogcneously-sclccte d samples (e.g., includin g 



radio-selected galaxies; [Alexander et al.l l2005bf ) . 



yielding results that conflict with X-ray studies 
of purely submil li meter-selected SMG sample s 
( Laird et al. 20101 : Georgantopoulos et al. 2011 ). 
To disentangle the relationship between SMGs 
and X-ray selected AGN, we need to overcome 
the uncertainty introduced by inhomogeneously 
selected samples, requiring X-ray spectral analyses 
of large, flux- limited samples of (sub) millimeter- 
selected SMGs with robust counterparts. 

In this work, we report on an X-ray stacking 
analysis of a sample of 38 SMGs detected in a 
1.2 mm map of the Lockman Hole North (LHN), 
one of the fields in the Spitzer Wid e- Area Infrared 
Extragalactic (SWIRE) Survey ( Lonsdale et al.l 



2003), using data from the Chandra-SWIRE sur 



vev (|Polletta et al.ll2006HWilkes et al.ll2009[) . The 
high radio counterpart iden tification rate of the 
LHN SMG sample (93%; iLindner et all l201l[ ) 
is afforded by t he extremely d eep 20 cm map of 
the same field (|Owen fc Morrison 2009.). and al- 
lows for reliable X-ray phot ometry. The sample 
benefits from spectroscopic ( PoUetta et al.l 20061 
lOwen fc Morrison! 120091: iFiolet et al.l l2010l) and 
optic ally-derived photometric (jStrazzuUo et al 



2010l ) redshifts. Additionally, analyses of Hersche l 



observations of the LHN ( Magdis et al. 2010t 



Roseboom et al ] l2012h have delivered reliable pho- 



tometric redshifts and infrared luminosities for a 
large fraction of the sample by fitting far-IR pho- 
tometry with thermal-dust spectral energy distri- 
bution (SED) models. 

In §2, we describe the observations used in our 
analysis. §3 outlines our X-ray stacking technique, 
and our method for deriving rest-frame luminosi- 
ties. In §4, we compare our results to previous 
X-ray studies of SMGs, and discuss the possible 
origins of the FeKa emission seen in our stacked 
spectrum. In §5, we present our conclusions. In 
our calculations, we assume a WMAP cosmol- 
ogy with i Jn = 70kms~^Mpc~\ = 0.27, and 
= 0.73 (jKomatsu et al.ll201ll ). 



2. Data and Sample Selection 

2.1. Millimeter observations 
and stacking sample 

Our SMG sample consists of 38 o f the 41 signif 



leant detections in the 1.2 mm map (jLindner et al 



20111 ) of the LHN made using the Max Planck 



Millim eter Bolometer (MAMBO; iKrevsa et al 
19981 ) array on the Institut de Radioastronomie 
Millimetrique 30 m telescope. We exclude one 
source that lacks a plausible 20 cm radio counter- 
part (L20), one that has a likely X-ray counter- 
part (L26), and one nearby galaxy at z = 0.044 
(L29) from the stacking sample. Of our final sam- 
ple of 38 galaxies, 37 (97%) have robust 20 cm ra- 
dio counter parts with a chance of spurious asso- 
ciation (P; iDownes et allll986l ) of P < 0.05; the 
remaining galaxy, L32, has P = 0.056. Stacking 
is performed with the coordinates of the SMGs' 
radio counterparts, which have a mean offset of 
2.4" with respect to the SMG centroids. Five of 
our stacking targets have posi tions that are not 
listed in the 20 cm catalog of I Owen fc Morrison 
()2008l) because they had S/N< 5.0 (L9, L28, 
and L36), or they were blended together with 
nearby r adio sources (L17 and L 3 9) during ex- 
traction (|Owen fc MorrisonI l2008t iLindner et al 



20111 ). The sample has a mean redshift of {z) = 2.6 
(see Table [1]). 

2.2. Chandra ACIS-I Observations 

Our X-ray data are from the 3 x 3-pointing 
raster mosaic of the LHN obtained with the 
Ad vanced CCD Imagin g Spectrometer (ACIS- 



I; IWeisskopf et"anil996l) o n the Chandra X-ray 



telescope bv iPoUetta et al.l (|2006l ). The final mo- 



saic comprises nine 70 ks pointings arranged with 
^ 2' overlap (sec Figure [1]). It covers a total 
area of ~ 0.7 deg^ and has a limiting conven- 
tional broad band [Be; 0. 5-8.0 keV) sensitivit y 



of ~ 4 X lO-^^erg s-^cm-^ (jPoUetta et al.ll2006[ ) 



Fiolet et al.l (|2009[ ) used these same data to search 
for a stacked X-ray signal among 33 Spitzer 
24 /im-selected starburst galaxies, and found no 
significant 0.3-8 keV-band emission. 

Within the sample of 41 MAMBO detec- 
tions in the LHN, only L26 has a likely X- 
ray counterpart (CXOSWJ1 4523.6 -h585601) in 
the catalog of IWilkes et al] (|2009l ). This X- 
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ray source has conventional broad band {Be; 
0.5-8.0 keV), soft band (Sc; 0.5-2.0 keV), and 
hard band {He; 2.0-8.0 keV) X-ray fluxes of 
fsc = 2.53 X 10-15 erg s^^cm-^ fsc = 1-21 x 
10~^^ ergs~^ cm~^ , and f He = 1.57x lO^^^ ergs^^ cm" 
respectively, and a hardness ratio of HR = 
—0.32^0 34. The hardness ratio is defined by 
HR = {He - Sc)/{Hc + Sc), where He and 
Sc are the counts in the Chandra conventional 
hard and soft bands, respectively. 

3. Stacking Analysis 

In this section we describe our reduction of the 
Chandra X-ray data products and the methods 
used in our stacking analysis. We use two tech- 
niques: (1) image-based stacking in binned X- 
ray maps f ^3.ip . and (2) a photon-based spec- 
tral stacking procedure using optimized apertures 
( §3.2p . The following subsections describe our im- 
plementation of these two methods. 

3.1. Image-based stacking 

We generate 1" x 1 "-pixel gridded maps of 
the total counts and effective exposure time in 
the Be, Se, and He energy bands using the 
Chandra Interactive Analys is of Observations 
(CIAO; iFruscione et al. 2006 ) script fluximage. 
The characteristic energies input to fluximage to 
compute effective areas were 4.00 keV, 1.25 kcV, 
and 5.00 kcV for the Be, Se, and He bands, 
respectively. Wc then used the CIAO scripts 
reproject_image to merge the maps of each ob- 
servation into one mosaic, and dmimgcalc to pro- 
duce an exposure-corrected flux image in units of 
[photons cm^^s^i]. 

Figure [2] shows the resulting stacked image in 
each energy band. The 40" x 40" S/N postage 
stamp images are shown with a color stretch from 
S/N=-4 to +4. The peak S/N is 3.2, 4.8, and 
2.0 in the Be, Se, and He bands, respectively. 
The peak of the strong stacked detection in the 
Se band has an offset from the mean radio coun- 
terpart ccntroid position of < 1". 

3.2. Optimized broad-band stacking 

Our second stacking technique does not use 
a binned X-ray map. Instead, we compute the 
stacked count rate and flux by directly counting 



photons at the stacking positions. The photomet- 
ric aperture at each stacking location is derived 
using a technique simil ar to the optimized sta cking 



sup- 



algorithm presented in lTreister et al.l ()2011 
plementary information) . 

The size and shape of the aperture at each 
stacking position is chosen to maximize the point- 
source S/N at that position on the ACIS-I chips. 
The apertures are constructed as follows. For each 
stacking position (shown in Figure [TJ, we (1) use 
the CIAO script mkpsf to generate a 2D image 
of the local Chandra PSF, (2) convolve this PSF 
with a Gaussian smoothing kernel (see below), and 
(3) find the enclosed-energy fraction (EEF) con- 
tour Ceef that maximizes the S/N of the fiux 
within the aperture. The area enclosed by this 
contour defines the aperture. Because Poisson 
noise from the X-ray background is stronger than 
the flux at each position and the average expo- 
sure time does not change rapidly with increasing 
aperture size, the S/N within the aperture can be 
parametrized by Ceef as 



S/Noc 



EEF 



y/ ^(Ceef) 



(1) 



where A(Ceef) is the total enclosed by the 
contour Ce ef ■ This expre s sion is the same as that 
derived by iTreister et al. ( 2011 ) , except that in- 
stead of using only circular apertures, we allow for 
non-axisymmetric apertures that follow the local 
shape of the Chandra PSF. 

We compensate for the change in shape of the 
Chandra PSF with photon energy by generating 
two optimal apertures at each stacking position, 
one for each energy band (i.e., for characteristic 
energies of 1.25 keV and 5 keV). Parts of the LHN 
Chandra mosaic were imaged multiple times due 
to the overlapping edges of the individual expo- 
sures (see Figure [1]). For these positions we flnd 
the total effective aperture by maximizing Equa- 
tion [1] using a linear combination of PSFs, one for 
each overlapping observation. 

We broaden the local Chandra PSFs to ac- 
commodate photons that do not lie at the stack- 
ing centers due to intrinsic wavelength offsets in 
the galaxies and astrometric errors. Previous X- 
ray stacking analyses find t he optimal aperture 
radius to be a = 1.25 --3.0" (|LeImiereLalJ[200^ 
Georgantopoulos et al.l 120111 ) . We find similarly 
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that our broad band stacked S/N is maximized 
with a smoothing kernel radius (Tkornci = 1-3" (see 
Figure [3]) , so wc adopt this value for our subse- 
quent broadband photometry. The smallest an- 
gular separation of any pair of stacking targets 
(15" for L17 and L39) is larger than the maximum 
radial extent of the largest X-ray stacking aper- 
ture, so we can ignore the effects of X-ray blending 
within our sample. 

The photons used for background subtraction 
are collected from arrays of large circular apertures 
positioned next to each stacking target position. 
The apertures are manually positioned to exclude 
any bright nearby X-ray sources that could con- 
taminate the background estimate. To avoid pos- 
sible systematic uncertainties assoc iated with the 



background subtraction (see, e.g., iTreister et al.l 
[20 11; .Willott .,20 ill) , we do not impose any S/N- 
based clipping or additional filtering in the back- 
ground regions. 

Table |3] shows the average stacked count rate 
and energy flux per galaxy in the three broad en- 
ergy bands. We find a significant stacked detec- 
tion in the soft band, and no significant detec- 
tion in the hard band. To convert the stacked 
count rate into energy flux, we used the web-based 
CIAO Portable Interactive Multi-Mission Simula- 
toiQ (PIMMS- version 4.4; Cycle 5). The fluxes 
are corrected for Galactic absorption using the col- 
umn density in the dire ction of the LHN center, 
TVh = 6.6 X 10^9 cm-2 (jStark et al.lll992[) . Our 
non-detection in the hard band leaves our calcula- 
tion of the hardness ratio relatively unconstrained, 
HR = — 0.68lo'32 (setting a limit on the photon in- 
dex r > 1.2), although it is clear that our sample 
has a steeply declining photon spectr um character- 



istic of star-forming galaxies (e.g., iRanalli et al 
2OI2I) . This estimate of HR was made after sub- 



tracting out the count rate in the soft band that 
is due to the strong FeKa line (see WS.'d^ . which 
is a ^ 20% contribution for 1.3"-broadened pho- 
tometric apertures. 

A high ph oton index o f T = 1 .6 (HR ^ -0.37) 
is found by iLaird et~aL (|2010l) . who stack on 
SCUBA-detected SMGs in the Chandra Deep 
Field North (CDF-N). An even steeper photon 
ind ex of F = 1.9 (HR :^ —0.4 9) is measured 
by ICeorgantopoulos et"~al (|20ll . who stack on 



LABOCA-detected SMGs in the Extended Chan- 
dra Deep Field South (ECDF-S). We have used 
values of F = 1.6 and F ~ 1.9 to compute the 
stacked flux in each energy band (e.g., see Table 
[3]), although the difference between the two es- 
timates is less than the Poisson uncertainty (see 
Table El). 

3.3. Optimized spectral stacking 

In addition to stacked broad-band fluxes, we 
also calculate the observed-frame and rest-frame 
stacked count-rate spectra for our sample. 

For each stacking target, we use redshift infor- 
mation in the following order of p riority, subject 
to avail ablility: (1) s pectroscopic J PoUetta et al 



i2006; Owen fc Morrison 2009: Fi olet et al.l 
(2) gersc/te/-based photome tric ( Magdis et al.l 



2010t iRoseboom et al.l |2012|). (3) A A-gualitv 



optical-based photometric ( Strazzullo et al. 201(J), 
and (4) m illimeter /radio photom etric estimated 



using the Carilli fc Yun ( ~999l) spec tral index 



„,20 cm 
"850 fim 



technique (jLindner et al. 2011 ). The red- 
shift distribution of our sample is shown in Figure 

SI 

We use a flat sum of the observed counts at each 
stacking-target position with no weighting factors. 
Although this technique gives more weight to the 
brightest members of the stack, it is necessary 
given that none of our stacking targets are indi- 
vidually detect ed and therefore S/N -based weights 
(used in, e.g., ITreister et al.l 1201 II ) cannot be re- 
liably assigned. For the rest-frame data, wc sep- 
arately coadd, blueshift, and bin the background 
photons to avoi d creating art ificial spectral fea- 



tures (see, e.g., lYaooob II2006I ). 



The uncertainty in the rest-frame energy of the 
photons Ai^rest as a function of the observed pho- 
ton energies i?obs due to the typical redshift error 
A2 is estimated by the equation 



AEr, 



E. 



obs 



(2) 



^http: / / cxc.harvard.edu/toolkit/pimms.jsp 



This uncertainty is always larger than the en- 
ergy resolution of the ACIS-I chips, so we set 
the rest-frame energy bin widths to match Ai?rest 
(Equation [2]) using our sample's average redshift 
(z) = 2.6 and redshift uncertainty (Az) = 0.4 (see 
Table E]). 

Figures [5] and |6| show the net observed and 
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rest-frame count-rate spectra for our SMG sam- 
ple, respectively. The rest-frame spectrum con- 
tains a 4cr emission feature with a centroid near 
6.7 keV, which we attribute to FeKa line emission 
from a mixture of Fe ionization states including 
FeXXV (see ii5.1.3p . It is apparent from this rest- 
frame spectrum that a significant fraction of the 
observed soft-band flux is due to the FeKa emis- 
sion line. If strong unresolved FeKa emission is 
a common feature in the X-ray spectra of other 
SMG samples, it may artificially lower their mea- 
sured HR values by inflating their observed soft- 
band fluxes. 

To ensure our stacking signal is not the result of 
contamination from a few strong targets, we per- 
formed a bootstrapping Monte Carlo analysis to 
recover the probability distribution for the single 
6.7keV-encrgy bin (bing.r). The mean number of 
on-target counts in bing.y is 16, while the mean 
number of background counts in the bin is 8. Fig- 
ure [7] shows that the resulting distribution closely 
matches that of an ideal Poisson distribution with 
a mean of 16, confirming that our stacking signal 
is characteristic of the entire sample, not a few 
outliers. 

3.4. Estimating Lx and LpcKa 

The mean stacked rest-frame hard-band X-ray 
luminosity, , of our sample is given by 

(Lhc) = ^TT fHo (dl) , (3) 

where Jhc is the stacked energy flux per galaxy in- 
side the observed 0.56-2.22 keV energy band (the 
2.0-8 keV energy band, redshifted by (z) = 2.6), 
and c?L is the luminosity distance at the redshift 
of each stacking target. We convert the observed 
count rate to an energy flux using PIMMS. The 
resulting rest frame X-ray luminosity is {Lhc) = 
(3.0 ± 1.1) X 1042gj.gs-i^ 

To estimate the equivalent width and line flux 
of the Fe Ka feature, we assume that the line emis- 
sion is contained only within the single elevated 
bin at 6.7 keV (see Figure[6]) and estimate the local 
count-rate continuum around the Fe Ka feature by 
averaging together the 8 bins between 3-9 keV (ex- 
cluding the bin containing the line). This results 
in an equivalent width of EW = 3.9 ± 2.5 keV. 
Although the EW is relatively unconstrained, it 
is > 1 kcV with 90% confldence. Using the nom- 



inal equivalent width and Equation IA7[ we find 
a mean stacked FeKa line flux of (/fcKq) — 
2.1 X 10~^^ergcm~^ s~^, and a mean line lumi- 
nosity (LpcKa) = (1.3 ±0.4) X lO^^g^gg-l^ 

Figure [3] also shows that that the S/N of the 
FeKa signal only drops at a larger radius than 
the broad-band signal. If we interpret the X-ray 
continuum as originating from the galaxies' nu- 
clear regions, then this relative offset between the 
Fe Ka emission and the X-ray continuum indicates 
that the Fe Ka photons in our sample are system- 
atically offset from the galaxies' centers. By mea- 
suring the distance between the peaks of the two 
curves in Figure [3l we estimate the radial offset 
to be ~ 1". For our computation of the FcKa 
line luminosity, we adopt an aperture broadening 
kernel suited to maximize the S/N of the FeKa 
emission line, (ikcrnci ~ 2.4". 

We use a two-sample Kolmogorov-Smirnov 
(KS) test to determine the signiflcance of the 
apparent angular extension of the FeKa emis- 
sion relative to the continuum emission. First, 
we compute the Chandra PSF at the location 
of each stacking target, then sample the PSFs 
at the positions of their respective collections of 
optimally-selected photons (see ^3.3p . The PSFs 
are peak-normalized and smoothed by an amount 
^smooth to reflect intrinsic wavelength offsets and 
astrometric errors in the photon positions. We 
then compare the cumulative distributions of the 
observed soft-band (0.5-2.0 keV) continuum pho- 
tons (excluding those in the rest-frame Fe Ka bin) 
and the rest-frame Fe Ka photons using the two- 
sample KS test to determine with what confidence 
1 — p (p is the KS test significance) we can rule 
out the null hypothesis that the two samples are 
drawn from a common distribution. When us- 
ing all 38 stacking positions, we find a maximum 
confidence of 1 — p = 0.71 at Csmooth = 0.5" (63 
continuum coimts and 15 FeKa counts). When 
we use only the stacking positions that have > 1 
FeKa photon, the maximum confidence occurs 
at the same value of (ismooth but has a reduced 
1 — p = 0.45 (29 continuum counts and 15 FeKa 
counts). Therefore, the extension in the FeKa 
emission relative to the continuum emission in- 
dicated by Figure [3] is a > Icr (71% confidence) 
effect, whose significance is limited primarily by 
the small number of FeKa photons. 
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4. Obscuration and Star Formation Rate 



Two galaxies in ou r sample have Lyp [8 



1000 iim \ estimated fromlMagdis et alj (|2010[ ). and 
12 from iRoseboom et al.] ~( 20121) . For the remain- 



ing galaxies without SED fits, we estimate Lir by 
scaling the SED from the nearby, bolometrically- 
star formation dominated ULIRG, Arp 220: 



L 



IR 



r Arp220 



Si, 



in terms of 



CfArp220 



1.2 mm X 



dL{zo) 



1 + z 
1 + zo 



l + ZQ 

1 + z 



(4) 



(5) 



^Arp220 ^ 3 ^ 1012 the observed 1.2mm flux 
density S'i.2mm! luminosity distance d^, target red- 
shift z, and Arp 220 rcdshift zq = 0.018. The Ljr 
values for all stacking targets are presented in Ta- 
ble m the average value of our whole sample is 



(Lir) = (2.4 ±0.2) X lO'^^ergs"^ 



(6) 



The mean 20 cm radio luminosity density L20 cm is 
calculated using our sam ple's redshift distribution 
and 20 cm flux densities ( Owen fc Morrisonll2008 : 



Lindner et al.ll2011^ 



{L 



20 cm, 



= (2.5 ±0.3) X lO^^ergs^^Hz"^ 



(7) 



m 
of 



We estimate the average star formation rate 
o ur sampl e using the scaling relations of 



. Kennicutt 
(I2OO3I) 



(Hks) 



m 



the 



infrared and 
at radio wavelengths, giving SFRir ~ 
(1100 ± 100)Moyr-i and SFR^adio ^ (1400 ± 
200) Mq yr"-'^, respectively. These values are 
consistent with each other, but greater than 
the estimat e using the X-ra y scaling relation 



from .Vattakunnel et all (120121 ). SFRx ^ (500 ± 



300) Mp, yr ^. All three scaling relations assume a 



Salpeteij (119551 ) initial mass function with limiting 
masses of 0.1 and 100 Mq. The SFR estimated 
using the X-ray luminosity may be low due to in- 
trinsic absorption. We can derive a lower limit on 
the average absorbing column in our sample by 
computing how much obscuration is required to 
reduce the value of SFRx from an intrinsic value 
consistent with SFRir and SFRiadio • In this case, 
would require Nyi ~ 2.3 x lO'^"^ cm^^ based 



we 



on our observed flux in the 0.55-2.77 keV band 
assuming F = 1.6 and using (z) = 2.6. If we use 
this argument to estimate the unabsorbed X-ray 
luminosity, we find {Lhc) — 9.2 x lO^^ergs"^. 

5. DISCUSSION 

5.1. Comparison to previous surveys 

In this section, we compare our results to those 
of previous X-ray analyses of SMG samples from 



the CDFN (Alexander et al. 2005b: Laird et al 



2010) and the (E)CDF-S (jGeorgantopoulos et al 
201 ih . 



5.1.1. Detection rate 

With only one s ignificant X-ra y coun terpart 



in the LHN, the iLindner et al.l (|201l[ ) SMG 



sample has an X-ray detection rate of 2^2%- 
Alexander et al.l (|2005bl ) find a high X-ray detec- 
tion rate of 85^20% among SMGs and submillimeter- 
targeted radio galaxies (whic h constitute 7 0% o f 
their sample) in the CDFN. iLaird et al.l (|2010h 
find a lower detection rate of 45 ± 8% using 
their purely submillimeter-selected sample de- 
rived from the inhomogenous ly covered SCUBA 
super map jBorvs et al.l boOSi) . The LESS sam- 
ple of lCeorgantopoulos et al.l (j201l[ ) is also purely 
submillimeter-selected and has an X-ray detection 
rate of 11^3%. However, unlike the SCUBA su- 
permap, the LESS survey is produced with a single 
observing mode, and with uniform coverage. 

We can place these four surveys in a common 
framework if we ask what fraction of SMGs in 
each survey have X-ray counterparts above the X- 
ray detection threshold in the L HN. In this case. 



we find 11 of 20 (55li|%) for lAlexander et al 



()2005bl) . of 35 (Ot^%) for iLaird et all (l2010l). 



and 1 1 of 126 (9^1%) for iGeorgantopoulos et al 



(|201ll ). The latter two are in agreement with our 
sample in the LHN. These results indicate that a 
lower X-ray detection rate may be more character- 
istic of strictly submillimeter-detected SMGs from 
surveys made with uniform coverage. 



5.1.2. LpcKa iiR vs. L 



20 cm 



Figure [8] shows our sample's average X-ray 
(corrected only for Galactic absorption), ra- 
dio, and IR luminosities com pared to those of 
other stacked SMG samples ( Laird et al. 2010t 
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GeorEantODOulos et alj 2011), individually X-ray- 


detected SMGs ( Alexander et 


al.ll2005b: Laird et al. 


2OIOI: Georffantopoulos et al. 


20 111), and nearby 



LIRG s and ULIRGs (sample drawn from llwasawa et al 
20091) . Where available, we use the Lir value 



from Table A2 of IPope et all (|2006f ) for the 



SMGs from the CDFN For the twelve SMGs in 
Alexander et al. J 2005lJ) that are not in the cat- 
alog of IPope et al] (|2006[ ). we scale Lfir — > ^ir 



using the average conversion factor / for the 
eight SMGs common between the two samples, 
/ = 1.42. The 870 ^m-detected SMGs from the 
(E)CDF-S are plotted with Lir = 10-1000 ^m 



The lo cal LIRGs and ULIRGs from llwasawa et al.l 
( 2009t ) also have their LpiR (40-400 ^m) estimates 
scaled by / = 1.42. We al so show the total sam- 
ple luminosity average for iLaird et~aL 

I (l2010h . in- 
cluding the contribution from their stacked SMGs 
that were not individually detected in the X-ray. 
The average properties of our stacking sample 
are in agreement with the total luminosities of 
Laird et all (|2010f) . 

Figures [8] and [9] also indicate the AGN classifi- 
cation of each galaxy. Galaxies whose mid-IR or 
X-ray spectral properties are consistent with emis- 
sion produced entirely by star formation arc plot- 
ted in red; while those req uiring the presence of an 
AGN are shown in blue. iGeorgantopoulos et al 



(|201l[ ) divide their sample by using a probabilis- 
tic approach; those galaxies requiring the pres- 
ence of a torus-dust component in their mid- 
IR SED according t o an F -test are categorized 
as AGN . Laird et al. ( 2010l ) and Alexander et al 



( 2005bl ) separate out the AGN based on the most 
favored mode l of th eir X-ray spectra according 
I1979I) statistic. The sample of 
20091) is divided based on hardness 



to the Cash 



Iwasawa et al 



ratio. 

The division between AGN and non-AGN sys- 
tems can be roughly determined based on the 
X-ray scaling relations of purely star-forming 
galaxies in the local univ erse (IRanalli et al.ll2003 : 
Vattakunnel et al. 20121 ). shown as solid black 
lines. The average pro perties of our s t acking sam- 
ple lie very near the iRanalli et al.l (j2003[ ) rela- 
tion. Considering the substantial intr insic scatter 



in the spectral classifications of the iLaird et al 
()2010l) sample, our stacking sample also proba- 
bly contains a substantial fraction of both star- 
formation-only and AGN- required systems. 



5.1.3. Fe Ka emission properties 

The FeKa photons in our stacking sample 
may be more spatially extended than the con- 
tinuum photons by ~ 1" (sec H3.4p . Extended 
and misalig ned Fe Kg emission h as been observed 
i n Arp220 (llwasawa et all l2005l) and NGC 1068 
(lYoung et al.ll200lh . We may also be blending 
together the emission from multiple components 
of merging systems of which only one component 
has strong Fe K g; emission (like, e.g., Arp 299; 



Ballo et all 120041 ). 

The bin width in our stacked rest-frame X-ray 
spectrum, which is set by the rcdshift uncertainties 
of our SMC sample, is larger than the rest energy 
separation between Fe Ka emission from neutral 
and highly-ionized iron (0.3 keV); therefore, it is 
difficult to determine the average Fe ionization 
fraction in our sample. Close inspection of the 
photons near the rest-frame Fc Ka line (see Fig- 
ure [T0|) reveals a range of values between 6.4 keV- 
7.2 keV, with a local maximum at 6.7 keV. Given 
that the Fc line photons are contributed fairly 
evenly by the 38 targets in our stacking sample, 
and have been assigned to their bins based on 
a wide variety of redshift estimation techniques 
(spectroscopic, optical-photometric, Herschel- 
photometric, millimeter/radio-photometric), they 
are unlikely all to be systematically biased high 
or low. Therefore, a significant fraction of the de- 
tected FeKa photons likely originate from the 
highly- ionized species of FeXXV or FeXXVI. 
However, the ^ 10% rest-frame uncertainty in 
the energy of each photon implies an uncertainty 
in the ccntroid of the line profile of CTcontroid = 
(FWHM)/(SNR) ~ 398 eV, insufficient to deter- 
mine the relative fractions of each ionization state 
with certainty. 

Strong emission (EW = 1.8 ± 0.9) from highly 
ionized Fe Ka has b een observed in the nearby 
ULIRG Arp 220 b y llwasawa et al l (l2005t ) using 
XMM - Newton, llwasawa et all (|2009t) also find 
strong 6.7 keV emission (EW = 0.9 ± 0.3 keV) 
from the stacked spectrum of nearby ULIRGs 
(including Arp 220) that have no evidence of 
AGN emission (termed X-ray-quict ULIRGs). 
I Alexander et al.l pOOSbl ) detected strong (EW ~ 
1 keV) Fe K a emission in the stacked SMC spec- 
trum of the six SMGs in their sample with 
Nh > 5 X 10^^ cm~^ , and find that the line cen- 
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troid is between 6.7 keV and 6.4 keV, indicating a 
substantial contribution from highly ionized gas. 

In Figure [11] we compare the relation between 
Lkq and Lir in our sample with those for other 
individual systems and stacked samples with mea- 
sured Fe Ka line luminosities and bolometrically 
dominant energy sources that are well under- 
stood. The dashed line represents a linear slope 
between Lkq and Lir and has been normalized 
to NGC 1068, a nearby prototypical Seyfert II 
LIRG. Red symbols represent systems that do not 
have significant observed AGN bolometric contri- 
butions, like SMGs and local X-ray quiet ULIRGs; 
the blue symbols represent systems that have sig- 
nificant bolometric AGN contributions. Figure [TT] 
shows that the relative Fe Ka/infrared luminos- 
ity fraction, LKa/-^iR, increases with increasing 
LiR. If the FeKa emission is due to AGN activ- 
ity, then this result may be in agreement with the 
observed trend that LIRGs/ULIRGs tend to be 
increa singly AGN-dom inated with increasing Lir 
(e.g., iTran et al.|[200l . 



5.2. Origin of the FeKa emission 

This section discusses three possible physical 
origins for the Fe Ka emission detected in our 
stacked SMG sample: supernovae, galactic-scale 
winds, and AGNs. Because a significant frac- 
tion of our sample's FeKa emission likely origi- 
nates from the highly-ionized species FeXXV (see, 
e.g.. Figure [To]) and because evidence for highly- 
ionized FeKa e mission from other (U)L IRGs ex- 
ists at both high lA lexan der et al.ll2005b[ ) and low 
(e.g., Ilwasawa et a l. 2005) redshifts, the following 
sections focus on the origin of this high-ionization 
component. 

5.2.1. Supernovae 

Here we consider if the observed FeKa fea- 
ture can be attributed to X-ray luminous su- 
pernovae. X-ray observations of the supernova 
SN 1986J in the nearby spiral galaxy NGC 831 re- 
veal strong hard-band emissi on and a sign ificant 
6.7 keV (EW < 500 cV) line ()Houck et al.l [T998'). 
Supernova 1986J decayed in the 2-10 keV band 
as ~ from 1991 to 1996. We wiU take a 
conservative approach and use only the luminos- 
ity information in this time interval for our cal- 
culation. Given the X-ray luminosity and decay 



rate of SN1986J (|Houck et al.lll998l ). and assum- 
ing the star formation rate of our sample of order 
SFR — lO'^M0yr~^ giving a supernova rate of 
10 SNyr~^, we would expect ^ 50 X-ray luminous 
supernovae to be visible at any given time. The 
combined supernova X-ray luminosity is therefore 
-^HcSNR — 10"*^ ergs~^. Considering the fact that 
prior to 1991 SN 1986J was probably still dimming 
at a rate close to oc this calculation is an un- 
derestimate. Therefore, supernovae like 1986J can 
satisfy the bolometric requirements for explaining 
the hard X-ray emission and the FeKa line that 
we see in our stacked SMG sample. 

However, if the supernovae associated with 
massive star formation are visible, then so must be 
high-mass X-ray binaries given the short time life- 
times of massive stars. These systems would dom- 
inate the hard X-ray emission from star-forming 
regions, and w ould severely dilute t he Fe Ka emis- 
sion (see, e.g., Iwasawa et al. 2009t ). We therefore 
rule out X-ray luminous supernovae and super- 
nova remnants as the source of the Fe K emission 
in our sample. 

5.2.2. Galactic- scale winds 



As discussed in Iwasawa et al. ( 2005h . who 



consider the 6.7 keV emission line in Arp220, a 
starburst-driven galactic-scale superwind of hot 
gas is energetically plausible as the source of the 
Fe Ka emission. Large outfiows could also explain 
why the the FeKa line emission appears more 
extended than the X-ray continuum emission in 
our stacking sample. To explore this scenario, we 
used the X-ra y spectral-fitting package XSPEC 
(|Arnaudlll99"6h to model an absorbed diffuse ther- 
mal X-ray (zphabs * mekal) spectrum and esti- 
mate the gas mctallicity required to produce the 
strong high-ionization Fe Ka emission detected in 
our SMG sample. We computed the model EW 
values using the spectral window 6.35-7.05 keV, 
the same energy width as the bin containing the 
Fe Ka emission in our stacked rest- frame spectrum 
(Figure [6|), which includes all FeKa ionization 
states. We fixed the gas te mperature to Arp 220 's 



p IM S 

Ilooi), 



best-fit value kT = 7.4 keV (jlwasawa et al. 
the gas density to n = lcm~^, the redshift to 
our sample's average (z) = 2.6, and the obscuring 
hydrogen column density to 2.3 x 10^"^ cni~^ (fj4|). 
Both the Fe Ka line luminosity and the contin- 
uum intensity vary linearly with Z, allowing us to 
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express the relation between EW and Z as 
6.67 



EW = 



1 



rkeV, 



(8) 



where Z' = 5.29 Zq. EW is approximately pro- 
portional to Z for Z <C Z' and approaehes the 
constant value 6.67 keV for Z Z'. Beeause of 
this non-linear behavior, an abundance of 0.94 Z^q 
can produce EW = IkeV (90% confidence lower- 
limit) while a significantly greater abundance Z ~ 
7.5 Zq is needed to explain our nominal value 
EW ~ 3.9 keV. If a significant amount of our 
rest-frame 2-10 keV luminosity is from X-ray bi- 
naries, incapable of generating the observed line 
emission, then the required metallicity would be 
even higher. While the lower limit on our mea- 
sured EW can be explained by thermal emission 
from a diffuse ionized plasma, especially consid- 
ering the extreme enrichment taking place in sys- 
tems like SMGs, generating an EW with a value 
close to our nominal measurement would require 
an unrealistic degree of high-z enrichment. 

5.2.3. AGN activity 

AGNs hidden behind large hydrogen column 
densities may be responsible for the observed 
FeKa emission in our sample. The FeKa emis- 
sion line is the signature spectral feature of 
the reprocessed (r eflected) spectrum of an AGN 



(|Matt et al.l |2000[) . As the ionizing luminos- 



ity increases, so does the ionization fraction of 
the gas, shifting the dominant emission feature 
from 6.4 keV (neutral and intermediate ioniza- 
tion states) to 6.7 keV (hehum-like Fe XXV) and 
6.9 keV (hydrogen-like Fe XXVI). 

Some insight into the properties of SMGs 
can be gained from reviewing the well-studied 
FeKa emission properties of AGN and nearby 
ULIRGs. Strong (EW ~ 1 keV) 6.7 keV emis- 
sion has been observed in systems that are 
bolom etrically AGN-dominated like IRAS 00182- 
7112 dNandra fc Iwasawal |2007|) and NGC 1068 



( Young et al.l 200 ll )^. as well as systems that ap- 
pear to be ene r getica lly AGN-free, l ike Arp220 



(llwasawa et al.l I2OO5I ) and IC694 (jBallo et al 



20041 ). However, direct evidence of a black hole 
accretion disk has been observed in Arp 220 by 
Downes fc Eckart I (|2007t ) with the detection of a 
compact (0.19" x 0.13") 1.3 mm continuum source 
in the center of the west nucleus torus. This source 



has a surface luminosity of ^ 5 x 10^** i© kpc~^, 
which is energetically incompatible with being 
powered by even the most extreme compact star- 
bursts known. Only an accretion disk can be 
responsible for heating the dust. Highly ionized 
Fe Ka emission h as also been observed in the AGN 
systems Mrk 273 (iBalest ra et al.ll20 0"5h. NGC 4945 



( Done et al.l 120031 ). and NGC 6240 (jBoller et al 



20031 ). along with a neutral FeKa component. 

The narrow 6.4 keV "cold" FeKa emission 
line is a ubiquitous feature in the spectra of 
optically-s elected active galaxies out to high red- 
shift (e.g., Corral et al.ll2008 : Iwasa , wa et al.ll2011 



Falocco et al 



2niih . 



However, 



Iwasawa et al 



(j201l[ ) also find evidence for highly ionized FeKa 
emission in two subsets of their X-ray selected 
AGN sample: Type I AGN with the highest Ed- 
dington ratios, and Type II AGN with the highest 
redshifts. The subsamples with highly ionized Ka 
emission show no evidence of a broad line FeKa 
feature; therefore, the highly- ionized FeKa pho- 
tons probably do not originate from the accretion 
disk, but from more distant and tenuous outfiow- 
ing gas. This scenario may also explain why the 
Fe Ka photons in our sample appear spatially 
extended with respect to the X-ray continuum 
photons. 

A significant caveat is that it remains difficult 
to reconcile the power source required to produce 
offsets as large as 1" in the photon distribution of 
our sample's stacked FeKa emission (relative to 
the nuclear continuum; see i j3.4[) given the sam- 
ple's low average X-ray luminosity. For exam- 
ple, we can calculate the maximum radial distance 
out to which low-density gas can remain highly 
photoionized by a single ionizing source by as- 
suming that our sample's stacked infrared lumi- 
nosity is produced by deeply-buried AGNs, i.e., 
iion — iiR = 2.4 X 10'*^ergs~^. Using the ion- 
ization parameter ^ = Lioy^/nR^ (logC ^ 2.8 is 
requir ed for a significant F eXXV ionization frac- 
tion: Kallman et al.l [2004 ) with n = Icm""^, we 
find i?max = \/LmJnS, — 2.0 kpc. At our sample's 
average redshift of {z) = 2.6, this corresponds to 
a typical angular offset of 0.24". Angular offsets 
larger than 0.24". like those tentatively indicated 
by our sample (see Figure [3]), can be explained 
by SMGs that host multiple distributed ionizing 
sources. In particular, the radio continuum emis- 
sion (defining our stacking positions) might be 



9 



more closely associated with the X-ray continuum 
than with Fe Ka line emission in a complex, multi- 
component system. These results highlight the 
importance of resolving the sizes and morpholo- 
gies of SMGs with high- resolution (sub)millimeter 
imaging (e.g., iTacconi et aLll2006l) . 

If the high- ionization FeKa emission is ulti- 
mately due to star-formation processes (shocked 
gas from SNe), and the SFR is traced by the iiR,, 
then we should expect a linear relation between 
Lkq and Xir. If the systems with the highest Ljr 
have an infrared contribution from obscured AGN 
that are not also emitting Fe Ka photons, then we 
would expect a slope that is even less than unity. 
However, Figure [TT] shows that Xkq is relatively 
much more dominant in SMGs and high-z AGN 
than in their lower-luminosity, lower-redshift ana- 
logues. This distinction indicates that highly ion- 
ized Fe Ka emission cannot be explained solely by 
star-formation processes and is more likely to be 
the result of AGN activity. 



tinuum spectral index F. 

We find a tentative indication (71% confidence) 
that our sample's stacked distribution of FeKa 
photons is more spatially extended than that of 
the X-ray continuum. If confirmed by future stud- 
ies, this result can help determine the physical ori- 
gin of the prominent Fe Ka emission in SMGs. 

We thank Jack Hughes, Dieter Lutz, and Alain 
Omont for useful discussions. We also thank the 
referee for useful comments that have improved 
the quality of this manuscript. This work has 
been supported by NSF grant AST-0708653, and 
has made use of the NASA/IPAC Extragalac- 
tic Database (NED) which is operated by the 
Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National 
Aeronautics and Space Administration. 

Facilities: IRAM:30m, CXO. 



6. Conclusions 

We analyze the X-ray properties of a complete 
sample of SMGs with radio counterparts from 
the LHN. This sample's X-ray detection rate of 
2l2% is consistent with those for other uniformly- 
mapped, submillimeter-detected samples, consid- 
ering the depth of our X-ray data. The X-ray un- 
detected SMGs show a strong stacked detection 
in the Sc band, and no significant detection in 
the He band, simil ar to results from SMG stack- 
ing in the CDF-N (iLaird et al.l l2010l) and CDF-S 
(jGeorgantopoulos et al. I l201l[) . 

We also use the available redshift information 
of our SMGs to compute the rest-frame, stacked 
count-rate spectrum of our sample. The rest- 
frame spectrum shows strong (EW > 1 keV) emis- 
sion from FeKa, possibly with contributions from 
FeXXV and FeXXVI. A comparison with other 
high-ionization Fe Ka-emitting systems from the 
literature indicates that accretion onto obscured 
AGNs is the likely explanation for the strong 
FcKa emission line. In our sample, the FeKa 
emission is responsible for ~ 20% of the observed 
soft-band X-ray fiux. Therefore, if strong Fe line 
emission is a common feature in other SMG sam- 
ples, it would significantly decrease the measured 
values of HR and lead to overestimates of the con- 
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A. Detailed descriptions of calculations 

A.l. Optimized spectral stacking method 

We begin by labeling all the photons within the optimized apertures (see Section [O]) of all of the iV = 38 
stacking targets with the index i, and those within the background regions for the N targets i'. Ei is the 
energy of the i*'' photon, and Ti is the total effective exposure time in the mosaic at the position of i. The 
notation i E j refers to all the photons that have energies located the j*'' energy bin. The stacked mean 
coimt rate in the j*'' energy bin, Rj, is then 

where Ki is the aperture correction for the optimal aperture of the energy band of the i*'* photon. The 
background mean count rate in the j*'' bin is 

where ce is the ratio of the areas of the background region of the stacking position of the i'*'* photon, and 
of the optimal aperture of that stacking position. It follows that the expected number of background counts 
in the j*'' bin, Nj, is 

Vj = (T,),^^. X i?;. (A3) 

Wc use the double-sided 68%-confidcncc upper and lower limits, ?7j'^^, and 77^°* ( Gehrels 19861) . to 
compute the la count rate deviations in the j*'' bin due to the background, 0-'^'/'°^: 

I hi/low I 

^^.hi/iow ^ - ^^^^ 

Therefore, the net count rate density per galaxy in the j*'' bin, Rj, is 

Rj = ^%;^ [B-^keV-M, (A5) 

where AEj is the width of the j*'' energy bin. We calculate the corresponding rest-frame spectrum Rj'^'^' by 
binning the photons according to their rest-frame energies, £'[°''' = Ei{l + Zi), where Zi is the redshift of the 
stacking target associated with the photon. 

A. 2. Fe Ka energy flux 

Wc use the CIAO script eff2evt to tabulate the local effective area Ai {Ai'), and quantum efficiency Qi 
{Qi'), for each photon in the on-target (background) apertures. 

The mean stacked on-target and background photon fluxes in the j*'' bin, Fj, and arc then 



and 

Q 



respectively. 



N ^. T,, A,, Q, K, c,. ' ^^^^ 
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Fig. 1.— SMG positions inside the LHN. The 
filled eircles mark the locations of the stacking tar- 
gets used in this work and are labeled according 
to the SMGs' ID numbers fi-om TablefH The loca- 
tions of the source without a reliable radio coun- 
terpart (L20), the nearby galaxy at z = 0.044 
(L29), and the source with a likely X-ray coun- 
terpart (L26) are also shown even though they are 
not used in our stacking analysis. The greyscale 
image shows the relative Chandra effective expo- 
sure time across the field. 




Fig. 2. — Stacked X-ray images showing the S/N 
in the Be, He, and Sc Chandra energy bands. 
The white cross hairs mark the stacking center. 

This 2-column preprint was prepared with the AAS lAT^X The color Stretch is S/N=[-4,-|-4] . 

maeros v5.2. 
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Fig. 3. — Stacked S/N of the observed Sc and 
rest-frame FeKa line fluxes as a function of the 
optimal aperture broadening kernel size CTkcmci- 
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Fig. 5. — Stacked observed-frame X-ray spectrum 
of our SMG sample. Above: Count rate spec- 
trum. The solid line and filled region represent 
the stacked on-target and background count rates. 
The error bars show the Poisson 68% double-sided 
confidence intervals due to the background. The 
dashed lines mark the extent of the Sc and He en- 
ergy bands. Below: S/N of each bin in the above 
spectrum. 
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Fig. 4. — Redshift distribution for our sample of 
38 SMGs. The filled histogram represents the red- 
shifts used in this work, with (z) = 2.6. The red 
dashed line represen ts the values as presented in 
Lindner et al.l (|2011l ) , which relied more heavily on 
spectral index-based photometric redshifts. 
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Fig. 6. — Stacked rcst-framc X-ray spectrum 
of our SMG sample. Above: Count rate spec- 
trum. The solid line and filled region represent 
the stacked on-target and background count rates. 
The error bars show the Poisson 68% double-sided 
confidence intervals due to the background only. 
The arrow marks the location of the rest-frame 
bin that is centered at 6.7 keV. Below: S/N of 
each bin in the above spectrum. 
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Fig. 7. — Bootstrap-Monte Carlo analysis of the 
number of FeKa counts in the 6.7keV-bin of 
the rest-frame stacked spectrum (Figure |6]). The 
histogram represents the total number of Monte 
Carlo trials returning the given number of total 
counts when randomly selecting 38 stacking posi- 
tions from the original 38, with replacement. 16 
counts lie in the 6.7 keV bin of the real data, and 
the mean Monte Carlo result is 15.9 counts. The 
solid curve shows an ideal Poisson distribution 
with mean ^ = 16. The shaded region represents 
the background level in this bin (8 counts). 
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Fig. 8. — Log X-ray luminosity vs. log infrared 
(8-1000 /im) luminosity for our SMG sample com- 
pared to other stacked and individually detected 
SMGs, LIRGs, and ULIRGs from the literature: 
SMGs in the CDF-N from Alexander et al. (2005b, 
triangles) and Laird et al. (2010, upside-down tri- 
angles) , SMGs in the ECDFS from Georgantopou- 
los et al. (2011, circles), and nearby LIRGs and 
ULIRGs from Iwasawa et al. (2009, squares). Red 
and blue symbols represent galaxies with "star 
formation-only" or "AGN-required" X-ray spec- 
tral classifications, respectively. The line repre- 
sents the relatio n for star-forming g alaxies in the 



local Universe ( Ranalli et al.l 120031 ) . The black 



hollow symbols represent the average luminosi- 
ties of all LIRGs/ULIRGs (Iwasawa et al. 2009, 
square), and all SMGs (including the stacked con- 
tribution) of the SCUBA supermap detections in 
the CDFN (Laird et al. 2010, upside-down trian- 
gle). 
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Fig. 9. — X-ray luminosity versus 20 cm spectral 
power. The symbols and colors are the same as 
in Figure [8] The over-plotted solid and dashed 
lines show the correlation between Lx and L20 cm 
along with la errors f or star-forming galaxies 
(IVattakunnel et al.|[201^ . 
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Fig. 10. — S/N of the FeKa emission as a func- 
tion of bin position and width. The red, green, 
and blue hnes use bin widths of 600 eV, 700 eV, 
and 800 eV (the bin width used in the rest frame 
spectrum in Figure [S] is 698 cV and is determined 
by the typical rcdshift uncertainty in our sam- 
ple). The black arrows mark the energies of 
Ka photons from Fcl (6.4keV), FeXXV (6.7keV) 
and FeXXVI (7.0 keV). The horizontal error bar 
shows the formal uncertainty in the line centroid 
Cccntroid = ±398 eV Centered on 6.7 keV. 
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Fig. 11. — FeKa line luminosity Lkq ver- 
sus LiR for our results (filled red circle), 
stacked SiVIGs in the CDF -N (empty red cir- 



cle; [Alexander et al ] l2005bl) . and stacked high-z 



Type II AGN (empty blue circle; llwasawa et al 



20111). The sq uares represent the LIRG NGC 1068 
(lYoung et al-T l200lh . the HvLIRG IRAS 00182- 



7112 (jNandra fc Iwasawall2007l). ar id the merger 
component IC 694 ( |Ballo et al.ll2004l ). The dashed 
line represents a linear relation that is normalized 
to NGC 1068. Blue (red) symbols denote sam- 
ples/systems that do (do not) show strong AGN 
bolometric signatures. 
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Table 1 
LHN MAMBO SMGs 



ID 


Name 




RA" 
J2000 


Doc" 
J2000 


z 


z typo*" 


z rcforoncc^ 


log(LiR/Lo)<i 


LiR reference^ 


LI 


MM J104700.1 + 


590109 


161 


75083 


59 


018778 




S 


P06 


13.3 


R12 


L2 


MM J104627.1 + 


590546 


161 


61192 


59 


095778 


4.09 


P-IR 


R12 


12.9 


L12 


L3 


MM J104631.4 + 


585056 


161 


63112 


58 


848889 


1.8 


CY 


Lll 


13.0 


L12 


L4 


MM J104607.4 + 


585413 


161 


53050 


58 


903889 


4.4 


CY 


Lll 


12.7 


L12 


L5 


MM J104725.2 + 


590339 


161 


85592 


59 


060139 


3.00 


P-O 


SIO 


13.1 


R12 


L6 


MM J104638.4 + 


585613 


161 


66112 


58 


936806 


2.03 


S 


FIO 


12.9 


R12 


L7 


MM J104700.1 + 


585439 


161 


75054 


58 


911500 


2.74 


P-IR 


R12 


12.7 


L12 


L8 


MM J104633.1 + 


585159 


161 


63779 


58 


866417 


2.07 


P-IR 


R12 


12.8 


L12 


L9 


MM J104704.9 + 


585008 


161 


77071 


58 


835750 


3.9 


CY 


Lll 


12.9 


L12 


LIO 


MM J104622.9 + 


585933 


161 


59604 


58 


993000 


3.01 


P-IR 


R12 


12.7 


L12 


Lll 


MM J104556.5 + 


585317 


161 


48721 


58 


888556 


1.95 


S 


FIO 


12.9 


R12 


L12 


MM J104448.0 + 


590036 


161 


19829 


59 


009972 


2.16 


P-O 


SIO 


12.7 


R12 


L13 


MM J104609.0 + 


585826 


161 


53646 


58 


974583 


0.32 


P-IR 


R12 


12.3 


L12 


L14 


MM J104636.1 + 


590749 


161 


64937 


59 


130139 


2.26 


P-IR 


R12 


12.9 


R12 


L15 


MM J104728.3 + 


585213 


161 


86654 


58 


870583 


2.78 


P-IR 


R12 


13.3 


R12 


L16 


MM J104720.9 + 


585151 


161 


83633 


58 


864722 


3.24 


P-IR 


R12 


12.8 


L12 


L17 


MM J104610.4 + 


590242 


161 


54383 


59 


045000 


2.65 


P-IR 


R12 


12.7 


L12 


LIS 


MM J104655.7 + 


585000 


161 


73021 


58 


834444 


2.15 


P-IR 


R12 


12.9 


L12 


L19 


MM J104502.1 + 


590404 


161 


25833 


59 


067583 


4.1 


CY 


Lll 


12.6 


L12 


L20 


MM J104617.0 + 


585444 


161 


57058 


58 


913694 


> 4.6 


CY 


Lll 






L21 


MM J104530.3 + 


590636 


161 


37575 


59 


110083 


2.52 


P-IR 


R12 


12.6 


L12 


L22 


MM J104603.8 + 


590448 


161 


51763 


59 


080389 


1.44 


P-IR 


R12 


12.3 


R12 


L23 


MM J104641.0 + 


585324 


161 


67129 


58 


890500 


3.6 


CY 


Lll 


12.6 


L12 


L24 


MM J104500.5 + 


590731 


161 


25279 


59 


126000 


3.24 


P-O 


SIO 


12.7 


R12 


L25 


MM J104540.3 + 


590347 


161 


41679 


59 


063333 


3.5 


CY 


Lll 


12.6 


L12 


L26 


MM J104522.8 + 


585558 


161 


34833 


58 


933611 


> 5.0 


CY 


Lll 






L27 


MM J104702.4 + 


585102 


161 


76000 


58 


850861 


1.62 


P-IR 


R12 


12.7 


L12 


L28 


MM J104620.9 + 


585434 


161 


58804 


58 


909722 


3.8 


CY 


Lll 


12.5 


L12 


L29 


MM J104556.1 + 


590914 


161 


48129 


59 


154528 


0.044 


S 


O09 


8.0 


R12 


L30 


MM J104510.3 + 


590408 


161 


29329 


59 


068639 


0.71 


P-IR 


R12 


12.5 


L12 


L31 


MM J104624.7 + 


585344 


161 


60379 


58 


896444 


2.90 


P-IR 


R12 


12.7 


R12 


L32 


MM J104539.6 + 


585419 


161 


41579 


58 


906917 


2.40 


P-IR 


MIO 


12.57 


MIO 


L33 


MM J104535.5 + 


585044 


161 


39608 


58 


847139 


2.63 


P-IR 


R12 


12.8 


L12 


L34 


MM J104453.7 + 


585838 


161 


22346 


58 


978194 


2.46 


P-IR 


R12 


12.6 


L12 


L35 


MM J104717.9 + 


585523 


161 


82546 


58 


923833 


2.14 


P-IR 


R12 


12.9 


L12 


L36 


MM J104608.1 + 


590744 


161 


53421 


59 


129444 


4.5 


CY 


Lll 


12.5 


L12 


L37 


MM J104610.8 + 


585242 


161 


54575 


58 


879139 


1.72 


P-IR 


MIO 


12.81 


MIO 


L38 


MM J104444.5 + 


590817 


161 


18704 


59 


138361 


3.6 


CY 


Lll 


12.5 


L12 


L39 


MM J104611.9 + 


590231 


161 


55012 


59 


042583 


2.59 


P-IR 


R12 


12.5 


L12 


L40 


MM J104658.7 + 


590633 


161 


74304 


59 


109306 


0.78 


P-IR 


R12 


12.5 


L12 


L41 


MM J104600.7 + 


585502 


161 


50129 


58 


917944 


1.49 


P-IR 


R12 


12.6 


L12 



"Position of 20 cm radio counterpart 

^Redshift typo: S — spe ctroscopic, P-IR, — H erschel-hascd photometric, P-O — optical-based photometric, CY — estimated using 
the spectral index c^ssQ^ni iCarilli & Yun|l 19991) . 

"Re ferences: P06 = | Polletta et al] 1200611. O09 ^lOwen fc MorrisonI II2009I'). FIO =IFiolet et al.l 12010) . MIO ^ IMagdis et al.l 12010) . 
SIO = IStrazzullo et al'l hoiOD . Lll = lLindner et al.l l201lh . R12 = IRoseboom et al.l 120121) . L12 = This Work 

■^8-1000 Aim 
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Table 2 
redshift uncertainties 



Rcdshift typo 



N galaxies 



References 



Spectroscopic 

Optical-based photometric 

Infrared-based photometric 

Spectral index-based estimate 
(Carilli fc Yun 1999) 



3 

23 



0.2 
0.4 

0.6 



Polletta et al. (2006 ') 
Fiolet et al. (2010) 
Owen fc Morrison (2009) 



Strazzullo et al. (2010 



Magdis et al. (2010) 
Roseboom et al. (2012) 



Lindn 



■ et al. (2011) 



Table 3 
X-RAY Stacking Results 



Band Energy Net Rate Flux (F = 1.6) Flux (F = 1.9) Lx (F = 1.6) 

keV 10"*^ s"^ 10"^'' ergcm"^ s"^ 10" erg cm"^ s"^ 10*'^ergs"^ 



Sc 


0.5-2.0 




^- ' -1.2 


4 9+^-^ 




He 


2.0-8.0 


1 n + 3.2 
-^■^-2.9 


< 9.7 


< 9.1 




Be 


0.5-8.0 
0.55-2.22 
0.55-2.78 


^■^-3.5 

8.oj:-i 


101^3:9 

5.711:? 


9 0+=* '' 
^■"-3.4 

5.7ti:? 


3.0 ± 1.1 
3.6 ± 1.3 



Note. — Unabsorbed fluxes calculated assuming the given photon index with Galactic absorption 
only. Hyd rogen column den sity taken as that of the central Chandra pointing, Nh — 6.6 X 
10^^ cm~^ jStark et al.lfT992D . Luminosity calculation uses P — 1.6. 
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